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Aerospace  Vehicles 


High  Speed  Systems 


Jm 


Power  and  Control 


Turbine  Engine 


CTO 


Rocket  Propulsion 


Motivation 


*  Biological  systems  not  necessarily  designed  for  optimal  flight 

*  Engineered  systems  don’t  have  requirements  related  to 
feeding,  care  for  young,  etc. 

*  Should  we  be  attempting  to  mimic  natural  systems,  knowing 
that  they  are  not  optimized  for  flight? 

*  What  would  a  biological  system  look  like  if  optimized  only  for 
flight? 

*  Can  we  use  engineering  design  and  optimization  to  create  a 
“flight-only  estimate”  of  the  biological  system? 
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Objective 


*  Understand  how  to  mechanize  multi-jointed  MAV 
wings  for  perching  and/or  flapping  applications 


*  Develop  an  energy-based  design  framework  for  the 
solution  of  combined  multi-physics,  multi-objective 

problems  Stiffened  area 


Cambering 
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Technical  Challenges 


*  Design  tool  for  multi-physics  analysis  and 
optimization  under  unsteady  aerodynamic  load  is 
not  well  established 

*  Identification  of  wing  morphology  requirements  is 
not  well  understood 

*  Performance  measures  such  as  energy  and 
efficiency  measures  for  unsteady  aerodynamic  flight 
are  not  well  defined 

*  Passive  shape  control  to  maximize  energy  efficiency 
is  not  well  exploited 
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Approach 


*  Student  1  (AFIT)  will  focus  on  the  distribution  of  skin 
material  to  meet  performance  objectives  after 
selecting  four  snap  shots  of  a  bird  wing 
configuration  during  perch 

*  Student  2  (UD)  will  extend  the  scope  of  the  research 
to  include  active  shape  control  (mechanism 
synthesis)  in  addition  to  skin  material  distribution 


Joints 

(Flapping  win; 


p - Fingers 

(Adv  structures) 

Wing  membrane 
^  (Advanced  materials) 

I,  Intrinsic  muscles 
■(Fluid-structure  interact) 


Swartz  &  Brauer,  Brown  Univ 


wing  hair  cell  sensors 
(Flight  &  flow  control) 
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Wing  Skin  Structure  Design 


*  Configuration  Selection 


Eagle  Owl  in  Loiter,  Dash,  and  Flare  Configurations 


Point  3 


Range 


Typical  Perching 
Trajectory  and 
Perching  Wing 
Configurations 
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Wing  Skin  Structure  Design 


*  Configuration  Selection 


Forward  Swept  Configuration  Zero  Sweep  Configuration 
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Drag  Coefficient 


Wing  Skin  Structure  Design 


A 


k 


•  Force  Estimation 

—  Forces  were  calculated  in  MATLAB  Vortex  Lattice  code  called 

Tornado 

—  Zero-lift ,  flat-plate  drag  coefficient  estimated  by  Tornado 
—  Drag  coefficient  related  to  angle  of  attack 
—  Force  on  each  panel  split  into  four  components  and  applied  to  the 


nodes 


3D  panels,  collocation  points  and  normals. 


0 


Anqle  of  Attack 


Viscous  Drog 


Body  x-coordinate 


Estimotion  Curve 


Example  of  Tornado  Vortex  Panels  Output 
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Wing  Skin  Structure  Design 


*  Perching  Data 


Wing  Configuration: 

Point  1:  Back  Swept 
Point  2:  Dive 
Point  3:  Zero  Sweep 
Point  4:  Forward  Swept 
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Wing  Skin  Structure  Design 


*  Force  Estimation 

-  Induced  drag  is  highest  for  Point  3,  Aerodynamic  Data  for  Birdwing 
not  Point  4,  and  lowest  at  Point  2  Alon9  Perchin9  Trajectory 

—  Side  forces  have  minimal 
influence  on  resulting  topologies 

—  Lift  highest  for  Point  3 

—  Axial  body  force  pushes  wing 
forward 

—  Most  bending  loads  about  10 
times  the  membrane  loads 

—  Viscous  drag  is  lowest  at  Point  4, 
even  though  the  Point  4  is  at  a 
high  angle  of  attack 


Point  1 

Point  2 

Point  3 

Point  4 

Vet. 

[m/s] 

10 

10.41 

10.11 

2.19 

AO  A 

[1 

6 

3.75 

10 

50 

Drag 

[N] 

0.0170 

0.0033 

0.0790 

0.0450 

Side 

[N] 

0.0001 

0.0075 

-0.0007 

-0.0087 

Lift 

[N] 

0.459 

0.112 

1.241 

0.190 

Fx 

“[NT 

-0.0304 

-0.0040 

-0.1371 

-0.1208 

F* 

[N] 

0.00010 

0.00749 

-0.00000 

-0.00871 

Fz 

[N] 

0.458 

0.112 

1.230 

0.101 

cL 

[-] 

0.220 

0.070 

0.512 

1.701 

cD 

[-] 

0.0085 

0.0022 

0.0328 

0.3958 

Cy 

[-] 

0.0029 

0.0051 

-0.0003 

-0.0757 

Re 

[-] 

90054 

137712 

91412 

19987 

Cdo 

H 

0.0101 

0.0082 

0.0101 

0.0113 

Swct 

M 

0.0081 

0.0444 

0.0775 

0.0785 

-Drjttf 

[N] 

0.2308 

0.1077 

0.4410 

0.0885 

Normal 

[N] 

0.0248 

0.0070 

0.0700 

0.0078 

Axial 

[N] 

0.2355 

0.1075 

0.4343 

0.0509 
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Wing  Skin  Structure  Design 


•  Results  -  Point  1 


Combined 
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Wing  Skin  Structure  Design 


•  Results  -  Point  2 


Combined 
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Wing  Skin  Structure  Design 


•  Results  -  Point  3 


(a)  Volf  =  0.2 


(b)  Volf  =  0.3 


(a)  Volf  =  0.2 


(b)  Volf  =  0.3 


(c)  Volf  =  0,4 


(d)  Volf  =  0,5 


Combined 
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Wing  Skin  Structure  Design 


•  Results  -  Point  4 


(a)  Volf  =  0.2 


(b)  Volf  =  0.3 


(a)  Volf  =  0.2 


(a)  Volf  =  0.2 


(c)  Volf  =  0.4  (d)  Volf  =  0.5 

Combined 
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Wing  Skin  Structure  Design 


•  Summary 

—  In  general,  structural  members  support  the  leading  edge 

—  Membrane  solutions  resemble  truss-like  structures,  and  bending 
solutions  resemble  beam-like  structures 

—  Membrane  solutions  clearly  dominate  the  combined  loading 

—  When  the  viscous  drag  distributed  over  the  surface  of  the  wing  is 
not  considered,  hybrid  solutions  occur 

—  Secondary  features  include  straight  battens  in  membrane 
structures,  and  branches  in  bending  structures 

—  Membrane  solution  must  support  out-of-plane  loading,  so  discrete 
“truss”  members  must  function  like  spars 

—  The  topology  constantly  changes  at  different  points  along  perching 
trajectory  so  we  need  an  active  mechanism  to  reconfigure  at 
different  loading  conditions  — ►  Wing  mechanism  design 
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Previous  Research 
(Multiple  Configurations) 


•  Generic  Surveillance  UAV  with  three  configurations 

•  Loiter  (configuration  0  =  reference) 

•  High  lift  (configuration  1) 

•  Climb  (configuration  2) 


-  Truss 

-  Frame 

Push  Active  Telescope 

-  Inactive  Telescope 

PushActuator 
►  PinAttachment 

•  Sliding  Attachment 
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Wing  Mechanism  Design 


*  Developing  design  tool  for  energy-based 
optimization  of  structure  topology 

*  Currently  includes... 

—  Geometry  Generator 
—  Pre-Processor 
—  Structural  Analysis 
—  Optimization  Routine 
—  Aerodynamic  Analysis  (in  progress) 

—  Post-Processor  (in  progress) 
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Wing  Mechanism  Design 


•  Geometry  Generator/Preprocessor 

—  Includes  a  GUI  for  ease  of  use 


—  Creates  a  parametrically  defined  wing  geometry 


Facilitates  future  optimization  routines  that  could  update  body 
geometry 
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Wing  Mechanism  Design 


*  Box  Substructure  Description 


16 

Nodes 


4 

Nodes 
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Wing  Mechanism  Design 


*  Structural  Analysis 

—  Implements  Standard  finite  element  approach 

—  Uses  a  condensed  frame  element  with  rotational 
springs  on  each  end 

—  Reduces  DoFs  thereby  decreasing  computational 


time  and  simplifies  programming 
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Wing  Mechanism  Design 


*  Optimization  Routine 

—  Globally  Convergent  Method  of  Moving  Asymptotes 

•  Developed  by  Svanberg 

•  One  of  the  most  used  methods  for  structural  optimization 


—  Problem  Formulation 

•  Minimize: 


Minimize: 


—  Shape  Error  and  Actuator  Usage  /» =  Wx  £ (u*r9et  -  u(p)i )2  -  W2 £ p) 


Subject  to: 

—  Static  Equilibrium 
—  Stroke  Limit 


Subject,  to: 

A,  =  KU  -  F  =  0 
=  Ei  -  £Lx  <  0 
!f  =  £  p*  ~  Nf  <  o 

i€/t 


»€T 


max 


<  o 


—  Attachment  Stiffness  f+v  =  £  +  E  p- 

t€/vl  i€/.2 

—  +/-  Volume  Fraction  f-v  =  ~Y  p>~  E  p+ <  o 

i€/d  i€//2 


i€4 

Static  Equilibrium 
Stroke  Limt 

Attachment  Placement  Limit 
Volume  Fraction  Limit 
V olume  Fraction  Limit 
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Wing  Mechanism  Design 


*  Aerodynamic  Analysis  (in  progress) 

—  Extracting  Aerodynamic  Influence  Coefficient  (AIC) 
matrix  from  Tornado  for  use  in  a  static  aeroelastic 
analysis 

—  Coupling  aerodynamic  loads  and  structural 
deformation 

—  Leveraging  the  aeroelastic  deformation,  it  is  assumed 
a  reduced  use  in  energy  design  may  be  found 

•  Post-  Processor 

—  Clearly  displays  the  results  from  the  design  tool 
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Research  Plans  for  Next  FY 


*  Key  energy  metrics  and  efficiency  measures  for 
optimal  multi-physics  designs 

*  Design  methodology  to  determine  passive  and 
active  shape  control  for  efficient  vehicle  flight 
performance 

*  Comparison  of  engineering  and  evolutionary  optimal 
solutions  for  similar  systems 


24 

Approved  for  public  release;  distribution  unlimited. 


Backup 
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Approach 


*  Utilize  design  optimization  techniques  for  efficient 
design  of  aeroelastic  reconfigurable  systems 
incorporating  distributed  actuation  and  compliance 

*  Develop  flight  energy  and  efficiency  measures  for 
topology  optimization 

*  Provide  understanding  of  a  systematic  design 
process  for  a  bio-mimetic  vehicle  design  problem 

S 

S 

s  ^ 

*  Select  “snapshots”  of  vehicle  ^ 
in  perching  maneuver  at 
different  times 

*  Optimize  based  on  multiple 
load  conditions 

*  Identify  suitable  objective 
functions  to  produce  “good” 
designs 
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Approach 


*  Student  2  (UD)  will  extend  the  scope  of  the  research 
to  include  mechanism  design  scheme  in  addition  to 
skin  material  distribution 


(Adv  structures) 


Joints 

(Flapping  wing) 


""  Wing  membrane 
(Advanced  materials) 


f  Intrinsic  muscles 
(Fluid-structure  interact) 


m 
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Wing  Skin  Structure  Design 


A 


k 


•  Optimality  Criteria  Method 

—  OC  method  is  a  bisection  method  based  on  the  fact  that  the 
material  volume  is  a  monotonically  decreasing  function  of  the 
Lagrange  multiplier 

—  Stationarity  point  is  achieved  when  volume  constraint  is  satisfied 
—  Update  scheme  given  by: 


p£+ 1  =  min 


such  that  the  volume  constraint  satisfies 

JV 


Xa^+i(  a)_  v  =  0 


—  OC  method  closely  related  to  fully  stressed  design,  where  all 

elements  have  same  strain  energy;  not  exactly  the  case,  because  of 
SIMP  model 
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Previous  Research 
(Flexible  Skin  Design) 


Target  shape  | 

Example  Target  Shapes 

^  Notional  Substructure 


P  E 

E:ju(Pe)= 

Two  Phase  Material  Solution 


•  Two-step  topology  optimization 
process 

-  Step  1 :  distribution  of  bulk 
material  properties 

-  Step  2:  distribution  of  multi¬ 
phase  material 

a) 


1.6979 

0.6230 

0 

d  = 

0.6230 

1.8880 

0 

xle3 

0 

0 

0.5066 

(2) 


Target  Reduced  Stiffness 
Matrix 


1.7179 

0.6076 

0 

QH  = 

0.6076 

1.9021 

0 

xle3 

0 

0 

0.5184 

Reduced  Stiffness  Matrix  from 
Homogenization  Routine 


Turning  Theory  Into  Application 
Reducing  Design  Time 
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Ad  Hoc  Solution 


Wing  Skin  Structure  Design 


*  Finite  Element  Derivation 

—  Membrane  Element 


—  Bending  Element 

—  Combined  Membrane/Bending  Element 


Superimposed  membrane  and  bending  plate  models  to  form  6-dof  model 


{/m} 
a  x  i 

l*m] 

8x8 

[0] 

8  x  12 

[0] 

20  x  4 

{<4n} 

8  x  i 

{fb} 

12  x  1 

►  = 

[0] 

12  x  8 

[kb) 

12  x  12 

i 

R}  ► 

12  x  1 

{/*} 
t  4  x  1  j 

[0] 

L  4  x  20 

[0] 

4x4. 

R*} 

w  4  x  1  j 

Fictitious  stiffness  matrix  added  for  "drilling''  degrees  of  freedom  to  avoid  singularities 


r  y 

A/,i 

1.0 

-0.5 

-0.5 

-0.5 

f  \ 

&zl 

a/22 

►  =  aEV 

-0.5 

1.0 

-0.5 

-0.5 

i 

0,2 

> 

A/,3 

-0.5 

-0.5 

1.0 

-0.5 

0,3 

A/,4 

J 

-0.5 

-0.5 

-0.5 

1.0 

0,4 
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Wing  Skin  Structure  Design 


A 
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*  Topology  Optimization 

—  Minimizing  compliance  equivalent  to  maximizing  stiffness 

—  Compliance  is  equivalent  to  the  strain  energy  of  a  deformed 
structure 

—  Volume  constraint  is  added  to  avoid  infinite  stiffness 
—  Nested  compliance  minimization  optimization  statement: 


min  c(p) 


P 


S.t.  {p}  {fl}  I  0)  0<  pmin  ^  Pc  ^  Pman  ^  1, .  .  .  ,  A 


where  the  compliance  c  is  defined  by 


c(p)  —  {F}7  {</},  where  {</}  solves:  {</}  =  {F} 
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Wing  Mechanism  Design 


*  Geometry  Generator/Preprocessor 

—  Generates  varying  degrees  of  mesh  connectivity  for 
the  initial  ground  structure  topology 


Degree  of  Connectivity  =  1  Degree  of  Connectivity  =  2 
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Wing  Skin  Structure  Design 


*  Solid  Isotropic  Material  with  Penalization  (SIMP) 

—  Penalizes  intermediate  thickness  values,  driving  thicknesses  towards 
a  discrete  solution 


—  Thicknesses  are  penalized  by  raising  the  element  thickness  to  a 
power  greater  than  1  in  the  constitutive  matrix:  r  v  0 
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